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Abstract 

These lectures are intended as an introduction to the technique of path integrals 
and their applications in physics. The audience is mainly first-year graduate students, 
and it is assumed that the reader has a good foundation in quantum mechanics. No 
prior exposure to path integrals is assumed, however. 

The path integral is a formulation of quantum mechanics equivalent to the standard 
formulations, offering a new way of looking at the subject which is, arguably, more 
intuitive than the usual approaches. Applications of path integrals are as vast as those 
of quantum mechanics itself, including the quantum mechanics of a single particle, 
statistical mechanics, condensed matter physics and quantum field theory. 

After an introduction including a very brief historical overview of the subject, we 
derive a path integral expression for the propagator in quantum mechanics, including 
the free particle and harmonic oscillator as examples. We then discuss a variety of 
applications, including path integrals in multiply-connected spaces, Euclidean path 
integrals and statistical mechanics, perturbation theory in quantum mechanics and in 
quantum field theory, and instantons via path integrals. 

For the most part, the emphasis is on explicit calculations in the familiar setting 
of quantum mechanics, with some discussion (often brief and schematic) of how these 
ideas can be applied to more complicated situations such as field theory. 



*Lectures given at Rencontres du Vietnam: Vlth Vietnam School of Physics, Vung Tau, Vietnam, 27 
December 1999 - 8 January 2000. 
t rbmack@lps . umont real . ca 
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1 Introduction 



1.1 Historical remarks 

We are all familiar with the standard formulations of quantum mechanics, developed more 
or less concurrently by Schroedinger, Heisenberg and others in the 1920s, and shown to be 
equivalent to one another soon thereafter. 

In 1933, Dirac made the observation that the action plays a central role in classical 
mechanics (he considered the Lagrangian formulation of classical mechanics to be more 
fundamental than the Hamiltonian one), but that it seemed to have no important role in 
quantum mechanics as it was known at the time. He speculated on how this situation might 
be rectified, and he arrived at the conclusion that (in more modern language) the propagator 
in quantum mechanics "corresponds to" expiS/h, where S is the classical action evaluated 
along the classical path. 

In 1948, Feynman developed Dirac's suggestion, and succeeded in deriving a third formu- 
lation of quantum mechanics, based on the fact that the propagator can be written as a sum 
over all possible paths (not just the classical one) between the initial and final points. Each 
path contributes expiS/h to the propagator. So while Dirac considered only the classical 
path, Feynman showed that all paths contribute: in a sense, the quantum particle takes all 
paths, and the amplitudes for each path add according to the usual quantum mechanical rule 
for combining amplitudes. Feynman's original paper,[| which essentially laid the foundation 
of the subject (and which was rejected by Physical Review!), is an all-time classic, and is 
highly recommended. (Dirac's original article is not bad, either.) 

1.2 Motivation 

What do we learn from path integrals? As far as I am aware, path integrals give us no 
dramatic new results in the quantum mechanics of a single particle. Indeed, most if not 
all calculations in quantum mechaincs which can be done by path integrals can be done 
with considerably greater ease using the standard formulations of quantum mechanics. (It is 
probably for this reason that path integrals are often left out of undergraduate-level quantum 
mechanics courses.) So why the fuss? 

As I will mention shortly, path integrals turn out to be considerably more useful in 
more complicated situations, such as field theory. But even if this were not the case, I 
believe that path integrals would be a very worthwhile contribution to our understanding of 
quantum mechanics. Firstly, they provide a physically extremely appealing and intuitive way 
of viewing quantum mechanics: anyone who can understand Young's double slit experiment 
in optics should be able to understand the underlying ideas behind path integrals. Secondly, 
the classical limit of quantum mechanics can be understood in a particularly clean way via 
path integrals. 

It is in quantum field theory, both relativistic and nonrelativistic, that path integrals 
(functional integrals is a more accurate term) play a much more important role, for several 

^References are not cited in the text, but a short list of books and articles which I have found interesting 
and useful is given at the end of this article. 
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reasons. They provide a relatively easy road to quantization and to expressions for Green's 
functions, which are closely related to amplitudes for physical processes such as scattering 
and decays of particles. The path integral treatment of gauge field theories (non-abelian 
ones, in particular) is very elegant: gauge fixing and ghosts appear quite effortlessly. Also, 
there are a whole host of nonperturbative phenomena such as solitons and instantons that 
are most easily viewed via path integrals. Furthermore, the close relation between statistical 
mechanics and quantum mechanics, or statistical field theory and quantum field theory, is 
plainly visible via path integrals. 

In these lectures, I will not have time to go into great detail into the many useful ap- 
plications of path integrals in quantum field theory. Rather than attempting to discuss a 
wide variety of applications in field theory and condensed matter physics, and in so doing 
having to skimp on the ABCs of the subject, I have chosen to spend perhaps more time 
and effort than absolutely necessary showing path integrals in action (pardon the pun) in 
quantum mechanics. The main emphasis will be on quantum mechanical problems which 
are not necessarily interesting and useful in and of themselves, but whose principal value is 
that they resemble the calculation of similar objects in the more complex setting of quantum 
field theory, where explicit calculations would be much harder. Thus I hope to illustrate the 
main points, and some technical complications and hangups which arise, in relatively famil- 
iar situations that should be regarded as toy models analogous to some interesting contexts 
in field theory. 

1.3 Outline 

The outline of the lectures is as follows. In the next section I will begin with an introduction 
to path integrals in quantum mechanics, including some explicit examples such as the free 
particle and the harmonic oscillator. In Section 3, I will give a "derivation" of classical 
mechanics from quantum mechanics. In Section 4, I will discuss some applications of path 
integrals that are perhaps not so well-known, but nonetheless very amusing, namely, the case 
where the configuration space is not simply connected. (In spite of the fancy terminology, 
no prior knowledge of high-powered mathematics such as topology is assumed.) Specifi- 
cally, I will apply the method to the Aharonov-Bohm effect, quantum statistics and anyons, 
and monopoles and charge quantization, where path integrals provide a beautifully intuitive 
approach. In Section 5, 1 will explain how one can approach statistical mechanics via path in- 
tegrals. Next, I will discuss perturbation theory in quantum mechanics, where the technique 
used is (to put it mildly) rather cumbersome, but nonetheless illustrative for applications in 
the remaining sections. In Section 7, I will discuss Green's functions (vacuum expectation 
values of time-ordered products) in quantum mechanics (where, to my knowledge, they are 
not particularly useful), and will construct the generating functional for these objects. This 
groundwork will be put to good use in the following section, where the generating functional 
for Green's functions in field theory (which are useful!) will be elucidated. In Section 9, I 
will discuss instantons in quantum mechanics, and will at least pay lip service to important 
applications in field theory. I will finish with a summary and a list of embarrassing omissions. 
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I will conclude with a few apologies. First, an educated reader might get the impression 
that the outline given above contains for the most part standard material. S/he is likely 
correct: the only original content to these lectures is the errors.0 

Second, I have made no great effort to give complete references (I know my limitations); 
at the end of this article I have listed some papers and books from which I have learned the 
subject. Some are books or articles wholly devoted to path integrals; the majority are books 
for which path integrals form only a small (but interesting!) part. The list is hopelessly 
incomplete; in particular, virtually any quantum field theory book from the last decade or 
so has a discussion of path integrals in it. 

Third, the subject of path integrals can be a rather delicate one for the mathematical 
purist. I am not one, and I have neither the interest nor the expertise to go into detail about 
whether or not the path integral exists, in a strict sense. My approach is rather pragmatic: 
it works, so let's use it! 



^Even this joke is borrowed from somewhere, though I can't think of where. 



2 Path Integrals in Quantum Mechanics 



2.1 General discussion 

Consider a particle moving in one dimension, the Hamiltonian being of the usual form: 

The fundamental question in the path integral (PI) formulation of quantum mechanics is: 
If the particle is at a position q at time t — 0, what is the probability amplitude that it will 
be at some other position q' at a later time t = T1 

It is easy to get a formal expression for this amplitude in the usual Schroedinger formu- 
lation of quantum mechanics. Let us introduce the eigenstates of the position operator g, 
which form a complete, orthonormal set: 

q k) = g k) , (g'l q) = <^(?' - j dq \q) {q\ = 1. 

(When there is the possibility of an ambiguity, operators will be written with a "hat"; 
otherwise the hat will be dropped.) Then the initial state is |^(0)) = \q). Letting the state 
evolve in time and projecting on the state k'), we get for the amplitude A, 

A = {q'\ ^{T)) ^ K{q', T; q, 0) = {q'\ e^^^ \q) . (1) 

(Except where noted otherwise, h will be set to 1.) This object, for obvious reasons, is known 
as the propagator from the initial spacetime point {q,0) to the final point {q',T). Clearly, 
the propagator is independent of the origin of time: K{q', T + t;q,t) = K{q', T; q, 0). 

We will derive an expression for this amplitude in the form of a summation (integral, 
really) over all possible paths between the initial and final points. In so doing, we derive the 
PI from quantum mechanics. Historically, Feynman came up with the PI differently, and 
showed its equivalence to the usual formulations of quantum mechanics. 

Let us separate the time evolution in the above amplitude into two smaller time evolu- 
tions, writing e'^^^ — e-^H{T-ti)^-imi _ rpj^g amplitude becomes 

A = {q'\ e-'"(T-'^)e-'"'' \q) . 
Inserting a factor 1 in the form of a sum over the position eigenstates gives 

A = {q'\ e-^(^-*^) J dq^ e"^^*^ \q) 



= ldqiK{q',T;q,,ti)K{qi,h;q,0). (2) 

This formula is none other than an expression of the quantum mechanical rule for combining 
amplitudes: if a process can occur a number of ways, the amplitudes for each of these ways 
add. A particle, in propagating from q to q', must be somewhere at an intermediate time ti, 
labelling that intermediate position qi, we compute the amplitude for propagation via the 
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point qi [this is the product of the two propagators in (Q)] and integrate over all possible 
intermediate positions. This result is reminiscent of Young's double slit experiment, where 
the amplitudes for passing through each of the two slits combine and interfere. We will look 
at the double-slit experiment in more detail when we discuss the Aharonov-Bohm effect in 
Section 4. 

We can repeat the division of the time interval T; let us divide it up into a large number 
of time intervals of duration 6 = T/N. Then we can write for the propagator 



A = {q'\ (e 



iH5 



N 



\q) = (9 |e 



iH5 



A' times 

We can again insert a complete set of states between each exponential, yielding 

'iHS 



A 



dqN-i \qN-i) {qN-i\e~ 



lN-2 \qN-2) {qN~ 



iH5 



dq2 \q2) {q2\ e 

dqi--- dqN-i {q'\ e~'^^ \qN-i) {Qn 
■■■(gi|e-^^^|g) 

dqi - ■ ■ dqN_iKq^^q^_^Kq^_^^q^_^ ■ 



dqi\qi) (gi|e-^^^|g) 
\qN-2) 



^1 e 



K K 



(3) 



where we have defined go = Qi In = <?'• (Note that these initial and final positions are not 
integrated over.) This expression says that the amphtude is the integral of the amplitude of 
all A"-legged paths, as illustrated in Figure 1. 



q=q, 





5 25 38 (N-l)S T t 

Figure 1: Amplitude as a sum over all A^- legged paths. 

Apart from mathematical details concerning the limit when N —>■ 00, this is clearly going 
to become a sum over all possible paths of the amplitude for each path: 

A = ^path; 
paths 



where 



= J dqi--- dqN-i, Apath 



K K - - - K K 

^^qN,qN^i'^^qN-i,qN-2 '■132,91 -"^ 91,90 • 



paths 
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Let us look at this last expression in detail. 

The propagator for one sub-interval is Kq^^^^q. = e"*^*^ We can expand the 

exponential, since 6 is small: 

= {q,+i\(l-tHS-^H'5' + ...^\q,) 

= ilj+illj) - {(lj+i\H \qj) + o{6^). (4) 
The first term is a delta function, which we can writeQ 

qj) = - qj) = I ^e^^^^''^--"^). (5) 

J Ztt 

In the second term of (^), we can insert a factor 1 in the form of an integral over momentum 
eigenstates between H and \qj); this gives 

= -^sf^(^ + V{q,^i)) e'-^fe--^^), (6) 
J 2tt \2m J 

using {q\p) = expipq. In the first line, we view the operator p as operating to the right, 
while V{q) operates to the left. 

The expression is asymmetric between qj and qj+i] the origin of this is our choice of 
putting the factor 1 to the right of H in the second term of (^. Had we put it to the left 
instead, we would have obtained V{qj) in (0). To not play favourites, we should choose some 
sort of average of these two. In what follows I will simply write V{qj) where g,- = ^{qj + qj+i). 
(The exact choice does not matter in the continuum limit, which we will take eventually; 
the above is a common choice.) Combining and @, the sub-interval propagator is 

= / ^e-^fe--^) (l - ^5 (f^ + V{q,)^ + o{6^)^ 

^e*^j('?i+i-9j)e-»<5-yfe:gj)(^l _^ o(5^)). (7) 

There are such factors in the amplitude. Combining them, and writing qj = [qj+i — qj)/ 5, 
we get 

^path = / li 17^ expiSY, {pjqj - H{pj, qj)), (8) 
j=o j=o 

where we have neglected a multiplicative factor of the form (1 -|- o((5^))^, which will tend 
toward one in the continuum limit. Then the propagator becomes 

K = J dqi--- dqN-iApi,th 

= y n y n ^ exp i5 ^ (p^g^ - H^p^, q^)). (9) 
i=i j=o i=o 



'Please do not confuse the delta function with the time interval, 5. 
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Note that there is one momentum integral for each interval (A^ total), while there is one 
position integral for each intermediate position (A^ — 1 total). 

If N —>■ oo, this approximates an integral over all functions p{t), q{t). We adopt the 
following notation: 



K = J Vp{t)Vq{t)expi dt {pq - H{p, q)) . 



(10) 



This result is known as the phase-space path integral. The integral is viewed as over all 
functions p{t) and over all functions q{t) where g(0) = q, q{T) = q'. But to actually perform 
an explicit calculation, (|1^) should be viewed as a shorthand notation for the more ponderous 
expression (P), in the limit N —>■ oo. 

If, as is often the case (and as we have assumed in deriving the above expression), the 
Hamiltonian is of the standard form, namely H = p^ /2m + V{q), we can actually carry out 
the momentum integrals in (^. We can rewrite this expression as 



K 



N-l N-l 

n dqjex^-i5 ^ V{qj] 

j=l j=0 



3=0 j=0 

The p integrals are all Gaussian, and they are uncoupled. One such integral is 



dp 
— ( 

271 



^i5(pq—p'^ /2m) 



^ iSrnq'^/2 
2^5 



(The careful reader may be worried about the convergence of this integral; if so, a factor 
exp —ep^ can be introduced and the limit e ^ taken at the end.) 
The propagator becomes 



K 



N-l N-l N-l 

n dqjex^-i5 ^ V{qj) J] 

j=l j=0 j=0 

/ ^ \N/2 .N-l N-l 

n dqjexpiS J2 | 

j=l j=0 



\27ri5) 



mq 



m .mq", 
exp 10 — - 

2-ni5 2 

2 \ 



V{q, 



The argument of the exponential is a discrete approximation of the action of a path passing 
through the points go = 1iQ.ii' ' ' 1 Qn-Ii In = <i' ■ As above, we can write this in the more 
compact form 



K 



Vq{t)e 



iS[q(t)] 



(12) 



This is our final result, and is known as the configuration space path integral. Again, (|T2| 
should be viewed as a notation for the more precise expression (pH]) , as A^ — ^ 00. 



2.2 Examples 

To solidify the notions above, let us consider a few explicit examples. As a first example, 
we will compute the free particle propagator first using ordinary quantum mechanics and 
then via the PI. We will then mention some generalizations which can be done in a similar 
manner. 
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2.2.1 Free particle 

Let us compute the propagator K{q', T; q, 0) for a free particle, described by the Hamiltonian 
H = /2m. The propagator can be computed straightforwardly using ordinary quantum 
mechanics. To this end, we write 

K = (g'|e-^^^|g) 

= {q'\e-^^^y'- f^\p){p\q) 



27r 

^e-^TpV2™(g'b)(p|g) 

/TT 

^^-iT{p2/2m)+i{q'-q)p^ ^^3) 

2% 



The integral is Gaussian; we obtain 



\2mTj ^ ' 

Let us now see how the same result can be attained using Pis. The configuration space 
PI (0) is 

= i^-^oo I M j y n exp [(gA. - g;v-i) + {qN-1 - qN-2) + ■ ■■ 

+ (^2 - qif + {qi - qof], 

where qo = q and q^ = q' are the initial and final points. The integrals are Gaussian, 
and can be evaluated exactly, although the fact that they are coupled complicates matters 
significantly. The result is 



liin ( ^ ] ^''^ gim(q'-g)2/27V5 

Af->oo V27rziV5/ 



But Nb is the total time interval T, resulting in 



K 



\2'KlT) 



in agreement with (|14|) . 

A couple of remarks are in order. First, we can write the argument of the exponential as 
T ■ \m{^{cl — q)/TY, which is just the action S[q(\ for a particle moving along the classical 
path (a straight line in this case) between the initial and final points. 
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Secondly, we can restore the factors of h if we want, by ensuring correct dimensions. The 
argument of the exponential is the action, so in order to make it a pure number we must 
divide by h; furthermore, the propagator has the dimension of the inner product of two 
position eigenstates, which is inverse length; in order that the coefficient have this dimension 
we must multiply by h~^^^. The final result is 



(15) 



This result typifies a couple of important features of calculations in this subject, which we 
will see repeatedly in these lectures. First, the propagator separates into two factors, one of 
which is the phase expiS[qc]/h. Second, calculations in the PI formalism are typically quite 
a bit more lengthy than using standard techniques of quantum mechanics. 



\2mhTj 



2.2.2 Harmonic oscillator 

As a second example of the computation of a PI, let us compute the propagator for the 
harmonic oscillator using this method. (In fact, we will not do the entire computation, but 
we will do enough to illustrate a trick or two which will be useful later on.) 

Let us start with the somewhat-formal version of the configuration-space PI, (|T2P- 

K{q',T;q,0) = J Vq{t)e'^^^'^'^l 

For the harmonic oscillator, 

/■^ /I 1 

S[(l{t)] = dt {-mq^ - -muj^q 



The paths over which the integral is to be performed go from g(0) = g to q{T) = q' . To do 
this PI, suppose we know the solution of the classical problem, qdt): 

qc + uj^qc = 0, qdO) = q, qdT) = q' . 

We can write q(t) = qdt) + y{t), and perform a change of variables in the PI to y{t), 
since integrating over all deviations from the classical path is equivalent to integrating over 
all possible paths. Since at each time q and y differ by a constant, the Jacobian of the 
transformation is 1. Furthermore, since qc obeys the correct boundary conditions, the paths 
y(t) over which we integrate go from y{0) = to y{T) = 0. The action for the path qc{t)+y{t) 
can be written as a power series in y: 



S[qc{t) + y{t)] = dt {^-mql - -mu^qc^j + (linear in y) + dt i^-mf - - 



2 1 

muj y 



The term linear in y vanishes by construction: q^ being the classical path, is that path for 
which the action is stationary! So we may write S[qc{t) + yit)] = S[qc(t)] + S[y(t)]. We 
substitute this into (|T^, yielding 



K{q', T- q, 0) = e'^t'^^W] / Vy{t)e''^^y'^'^\ (16) 
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As mentioned above, the paths y{t) over which we integrate go from y{0) = to y(T) = 0: 
the only appearance of the initial and final positions is in the classical path, i.e., in the 
classical action. Once again, the PI separates into two factors. The first is written in terms 
of the action of the classical path, and the second is a PI over deviations from this classical 
path. The second factor is independent of the initial and final points. 

This separation into a factor depending on the action of the classical path and a second 
one, a PI which is independent of the details of the classical path, is a recurring theme, 
and an important one. Indeed, it is often the first factor which contains most of the useful 
information contained in the propagator, and it can be deduced without even performing a 
PI. It can be said that much of the work in the game of path integrals consists in avoiding 
having to actually compute one! 

As for the evaluation of (0), a number of fairly standard techniques are available. One 
can calculate the PI directly in position space, as was done above for the harmonic oscillator 
(see Schulman, chap. 6). Alternatively, one can compute it in Fourier space (writing y{t) = 
J2k (^k sin(/c7rt/T) and integrating over the coefficients {afc}). This latter approach is outlined 
in Feynman and Hibbs, Section 3.11. The result is 



The classical action can be evaluated straightforwardly (note that this is not a PI problem, 
nor even a quantum mechanics problem!); the result is 



We close this section with two remarks. First, the PI for any quadratic action can be 
evaluated exactly, essentially since such a PI consists of Gaussian integrals; the general 
result is given in Schulman, Chapter 6. In Section 6, we will evaluate (to the same degree 
of completeness as the harmonic oscillator above) the PI for a forced harmonic oscillator, 
which will prove to be a very useful tool for computing a variety of quantities of physical 
interest . 

Second, the following fact is not difficult to prove, and will be used below (Section 4.2.). 
K{q', T; q, 0) (whether computed via Pis or not) is the amplitude to propagate from one 
point to another in a given time interval. But this is the response to the following question: 
If a particle is initially at position q, what is its wave function after the elapse of a time T? 
Thus, if we consider as a function of the final position and time, it is none other than 
the wave function for a particle with a specific initial condition. As such, the propagator 
satisfies the Schroedinger equation at its final point. 




(17) 




3 The Classical Limit: "Derivation" of the Principle 
of Least Action 



Since the example calculations performed above are somewhat dry and mathematical, it is 
worth backing up a bit and staring at the expression for the configuration space PI, (|T2|): 



This innocent-looking expression tells us something which is at first glance unbelievable, and 
at second glance really unbelievable. The first-glance observation is that a particle, in going 
from one position to another, takes all possible paths between these two positions. This 
is, if not actually unbelievable, at the very least least counter-intuitive, but we could argue 
away much of what makes us feel uneasy if we could convince ourselves that while all paths 
contribute, the classical path is the dominant one. 

However, the second-glance observation is not reassuring: if we compare the contribution 
of the classical path (whose action is >S'[gc]) with that of some other, arbitrarily wild, path 
(whose action is we find that the first is expiS[qc] while the second is expiS[qm]. 

They are both complex numbers of unit magnitude: each path taken in isolation is equally 
important. The classical path is no more important than any arbitrarily complicated path! 

How are we to reconcile this really unbelievable conclusion with the fact that a ball thrown 
in the air has a more-or-less parabolic motion? 

The key, not surprisingly, is in how different paths interfere with one another, and by 
considering the case where the rough scale of classical action of the problem is much bigger 
than the quantum of action, h, we will see the emergence of the Principle of Least Action. 

Consider two neighbouring paths q{t) and q'{t) which contribute to the PI (Figure 2). 
Let q'{t) = q{t) +ri{t), with ri{t) small. Then we can write the action as a functional Taylor 
expansion about the classical path:^ 



S[q'] = S[q + v] = S[q] + J dtv{t)f^ + o{v 



5q{t) 

The two paths contribute expz5'[g]//i and expiS'[g']/^ to the PI; the combined contribution 
is 

A^e^^W/^ (^l + expl|rftr/(t)^^, 

where we have neglected corrections of order rj^. We see that the difference in phase be- 
tween the two paths, which determines the interference between the two contributions, is 
h-^ !dtr^{t)5S[q]/5q{t). 

We see that the smaller the value of the larger the phase difference between two 
given paths. So even if the paths are very close together, so that the difference in actions 
is extremely small, for sufficiently small % the phase difference will still be large, and on 
average destructive interference occurs. 

"^The reader unfamiliar with manipulation of functionals need not despair; the only rule needed beyond 
standard calculus is the functional derivative: 5q{t) / 5q{t') = 5{t — t'), where the last 6 is the Dirac delta 
function. 
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q 




t 



Figure 2: Two neighbouring paths. 

However, this argument must be rethought for one exceptional path: that which extrem- 
izes the action, i.e., the classical path, qc{t). For this path, S[qc + 1]] = S[qc] + 0(77^). Thus 
the classical path and a very close neighbour will have actions which differ by much less than 
two randomly-chosen but equally close paths (Figure 3). This means that for fixed closeness 




t 



Figure 3: Paths near the classical path interfere constructively. 

of two paths (I leave it as an exercise to make this precise!) and for fixed h, paths near the 
classical path will on average interfere constructively (small phase difference) whereas for 
random paths the interference will be on average destructive. 

Thus heuristically, we conclude that if the problem is classical (action ^ h), the most 
important contribution to the PI comes from the region around the path which extremizes 
the PI. In other words, the particle's motion is governed by the principle that the action is 
stationary. This, of course, is none other than the Principle of Least Action from which the 
Euler- Lagrange equations of classical mechanics are derived. 



4 Topology and Path Integrals in Quantum Mechanics: 
Three Applications 

In path integrals, if the configuration space has holes in it such that two paths between 
the same initial and final point are not necessarily deformable into one another, interesting 
effects can arise. This property of the configuration space goes by the following catchy name: 
non-simply-connectedness. We will study three such situations: the Aharonov-Bohm effect, 
particle statistics, and magnetic monopoles and the quantization of electric charge. 

4.1 Aharonov-Bohm effect 

The Aharonov-Bohm effect is one of the most dramatic illustrations of a purely quantum 
effect: the influence of the electromagnetic potential on particle motion even if the particle is 
perfectly shielded from any electric or magnetic fields. While classically the effect of electric 
and magnetic fields can be understood purely in terms of the forces these fields create on 
particles, Aharonov and Bohm devised an ingenious thought-experiment (which has since 
been realized in the laboratory) showing that this is no longer true in quantum mechanics. 
Their effect is best illustrated by a refinement of Young's double-slit experiment, where 
particles passing through a barrier with two slits in it produce an interference pattern on a 
screen further downstream. Aharonov and Bohm proposed such an experiment performed 
with charged particles, with an added twist provided by a magnetic flux from which the 
particles are perfectly shielded passing between the two slits. If we perform the experiment 



Impenetrable 
shield 



Interference 
pattern shifts 



Figure 4: Aharonov-Bohm effect. Magnetic flux is conflned within the shaded area; particles 
are excluded from this area by a perfect shield. 

first with no magnetic fiux and then with a nonzero and arbitrary fiux passing through the 
shielded region, the interference pattern will change, in spite of the fact that the particles are 
perfectly shielded from the magnetic field and feel no electric or magnetic force whatsoever. 
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Classically we can say: no force, no effect. Not so in quantum mechanics. Pis provide a very 
attractive way of understanding this effect. 

Consider first two representative paths qi(t) and q2(t) (in two dimensions) passing 
through slits 1 and 2, respectively, and which arrive at the same spot on the screen (Figure 
5). Before turning on the magnetic field, let us suppose that the actions for these paths are 
S[qi\ and S'[q2]. Then the interference of the amplitudes is determined by 



The relative phase is = ('S'[q2] — S'[qi])//i. Thus these two paths interfere constructively 
if 012 = 2mT, destructively if 4>i2 = (2n + l)7r, and in general there is partial cancellation 
between the two contributions. 




Figure 5: Two representative paths contributing to the amplitude for a given point on the 
screen. 



How is this result affected if we add a magnetic field, B? We can describe this field by 
a vector potential, writing B = V x A. This affects the particle's motion by the following 
change in the Lagrangian: 

L(q, q) ^ L'(q, q) = L(q, q) - -v A(q). 

c 



Thus the action changes by 



.i/*v.A(q)^-l/*MO.A«,W). 



This integral is J dq^- A(q), the line integral of A along the path taken by the particle. So 
including the effect of the magnetic field, the action of the first path is 



^'[qi] = ^[qi] - - / c^q-A(q), 

C Jqi(t) 



qi(«) 

and similarly for the second path. 

Let us now look at the interference between the two paths, including the magnetic field. 



giS'[qi]/7i _^ gi5'[q2l/?i 



gi5'[qi]/a 1^-^ _,_ gi(5'[q2]-S'[qi])/nj 



(18) 
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where the new relative phase is 

0'i2 = <Pi2-^(f dq- A(q) - / rfq ■ A{q)] . (19) 

But the difference in hne integrals in (|19D is a contour integral: 

/ dq ■ A(q) - / dq- A(q) = i dq- A(q) = 

Jci2{t) Jcil{t) J 

$ being the ffux inside the closed loop bounded by the two paths. So we can write 



= '^^^ ^ Ye 



It is important to note that the change of relative phase due to the magnetic field is 
independent of the details of the two paths, as long as each passes through the corresponding 
slit. This means that the PI expression for the amplitude for the particle to reach a given 
point on the screen is affected by the magnetic field in a particularly clean way. Before the 
magnetic field is turned on, we may write A = Ai + A2, where 

A^= f Pqe^^W/^ 

Jslit 1 

and similarly for A2. Including the magnetic field, 

Jslit 1 ' 

where we have pulled the line integral out of the PI since it is the same for all paths passing 
through slit 1 arriving at the point on the screen under consideration. So the amplitude is 



' Ii d.^-J^/^<^ {Ai + e^^^^^^'^A 



2 

2 

The overall phase is irrelevant, and the interference pattern is influenced directly by the 
phase e^/he. If we vary this phase continuously (by varying the magnetic flux), we can 
detect a shift in the interference pattern. For example, if e^/he = vr, then a spot on the 
screen which formerly corresponded to constructive interference will now be destructive, and 
vice- versa. 

Since the interference is dependent only on the phase difference mod 2it, as we vary the 
flux we get a shift of the interference pattern which is periodic, repeating itself when e^/he 
changes by an integer times 27r. 
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4.2 Particle Statistics 

The path integral can be used to see that particles in three dimensions must obey either Fermi 
or Bose statistics, whereas particles in two dimensions can have intermediate (or fractional) 
statistics. Consider a system of two identical particles; suppose that there is a short range, 
infinitely strong repulsive force between the two. We might ask the following question: if at 
t = the particles are at qi and q2, what is the amplitude that the particles will be at q'^^ 
and q'g at some later time T? We will first examine this question in three dimensions, and 
then in two dimensions. 

4.2.1 Three dimensions 

According to the PI description of the problem, this amplitude is 



where we sum over all two-particle paths going from qi, q2 to q'^, q'g. 

However there is an important subtlety at play: if the particles are identical, then there 
are (in three dimensions!) two classes of paths (Figure 6). 



Even though the second path involves an exchange of particles, the final configuration is 
the same due to the indistinguishability of the particles. 

It is more economical to describe this situation in terms of the centre-of-mass position 
Q = (cti + q2)/2 and the relative position q = q2 — qi- The movement of the centre of mass 
is irrelevant, and we can concentrate on the relative coordinate q. We can also assume for 
simplicity that the final positions are the same as the initial ones. Then the two paths above 
correspond to the paths in relative position space depicted in Figure 7. 

The point is, of course, that the relative positions q and — q represent the same configu- 
ration: interchanging qi and q2 changes q — — q. 

We can elevate somewhat the tone of the discussion by introducing some amount of 
formalism. The configuration space for the relative position of two identical particles is not 

— {0},| as one would have naively thought, but (R^ — {0})/Z2. The division by the factor 

^Recall that we have supposed that the particles have an infinite, short-range repulsion; hence the sub- 
traction of the origin (which represents coincident points). 



,i5[qi(i),q2(i)] 





Figure 6: Two classes of paths. 



4 TOPOLOGY AND PATH INTEGRALS 



18 



Z2 indicates that opposite points in the space — {0}, namely any point q and the point 

diametrically opposite to it — q, are to be identified: they represent the same configuration. 
We must keep this in mind when we attempt to draw paths: the second path of Figure 7 is 
a dosed one. 

A topological space such as our configuration space can be characterized as simply con- 
nected or as non-simply connected according to whether all paths starting and finishing at 
the same point can or cannot be contracted into the trivial path (representing no relative 
motion of the particles). It is clear from Figure 7 that the first path can be deformed to the 
trivial path, while the second one cannot, so the configuration space is not simply connected. 
Clearly any path which does not correspond to an exchange of the particles (a "direct" path) 
is topologically trivial, while any "exchange" path is not; we can divide the space of paths 
on the configuration space into two topological classes (direct and exchange). 

Our configuration space is more precisely described as doubly-connected, since any two 
topologically nontrivial (exchange) paths taken one after the other result in a direct path, 
which is trivial. Thus the classes of paths form the elements of the group Z2 if we define 
the product of two paths to mean first one path followed by the other (a definition which 
extends readily to the product of classes). 

One final bit of mathematical nomenclature: our configuration space, as noted above, is 
(R^ — {0})/Z2, which is not simply connected. We define the (simply connected) covering 
space as the simply connected space which looks locally hke the original space. In our case, 
the covering space is just R^ — {0}. 

At this point you might well be wondering: what does this have to do with Pis? We 
can rewrite the PI expression for the amplitude as the following PI in the covering space 
R3 - {0}: 

direct exchange 

= A(q,T;q,0)+l(-q,T;q,0). (20) 

The notation A is used to indicate that these Pis are in the covering space, while A is a PI 
in the configuration space. The first term is the sum over all paths from q to q; the second 
is that for paths from q to — q. 




Figure 7: Paths in relative coordinate space. 
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Notice that each sub-path integral is a perfectly respectable PI in its own right: each 
would be a complete PI for the same dynamical problem but involving distinguishable parti- 
cles. Since the PI can be thought of as a technique for obtaining the propagator in quantum 
mechanics, and since (as was mentioned at the end of Section 2) the propagator is a solution 
of the Schroedinger equation, either of these sub-path integrals also satisfies it. It follows 
that we can generalize the amplitude A to the following expression, which still satisfies the 
Schroedinger equation: 

^(q, T; q, 0) ^ A^q, T; q, 0) = ^ e^^^ + e'^ ^ e'^^ 

direct exchange 

= A(q,T;q,0) + e*<^A(-q,T;q,0). (21) 

This generalization might appear to be ad hoc and ill-motivated, but we will see shortly that 
it is intimately related to particle statistics. 

There is a restriction on the added phase, 0. To see this, suppose that we no longer insist 
that the path be a closed one from q to q. Then (plD generalizes to 

A^iq!, T- q, 0) = A(q', T; q, 0) + e*'^A(-q', T; q, 0). (22) 

If we vary q' continuously to the point — q', we have 

^<^(-q', T; q, 0) = A(-q', T; q, 0) + e'^A{c^, T; q, 0). (23) 

But since the particles are identical, the new final configuration — q' is identical to old one 
q'. (p2D and ([23|) are expressions for the amplitude for the same physical process, and can 
differ at most by a phase: 

A^q',T; q,0) = e^M<^(-q',T;q,0). 
Combining these three equations, we see that 

A(q', T; q, 0) + e^*A(-q', T; q, 0) = e^" (^(-q', T; q, 0) + e^<^A(q', T; q, 0)) . 
Equating coefficients of the two terms, we have a = (p (up to a 2tt ambiguity), and 

e'"^ = 1. 

This equation has two physically distinct solutions: = and = tt. (Adding 2mi results 
in physically equivalent solutions.) 
If = 0, we obtain 

A(q, T; q, 0) = A(q, T; q, 0) + A{-q, T; q, 0), (24) 

the naive sum of the direct and exchange amplitudes, as is appropriate for Bose statistics. 
If, on the other hand, = vr, we obtain 

Aiq, T; q, 0) = A(q, T; q, 0) - A(-q, T; q, 0). (25) 

The direct and exchange amplitudes contribute with a relative minus sign. This case de- 
scribes Fermi statistics. 

In three dimensions, we see that the PI gives us an elegant way of seeing how these two 
types of quantum statistics arise. 
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4.2.2 Two dimensions 

We will now repeat the above analysis in two dimensions, and will see that the difference is 
significant. 

Consider a system of two identical particles in two dimensions, again adding a short- 
range, infinitely strong repulsion. Once again, we restrict ourselves to the centre of mass 
frame, since centre-of-mass motion is irrelevant to the present discussion. The amplitude 
that two particles starting at relative position q = {q^, Qy) will propagate to a final relative 
position q' = (g^, q'y) in time T is 

^(q',r;q,0) = ^ e'^[mitU2it)]^ 

paths 

the sum being over all paths from q to q' in the configuration space. 

Once again, the PI separates into distinct topological classes, but there are now an 
infinity of possible classes. To see this, consider the three paths depicted in Figure 8, where 
for simplicity we restrict to the case where the initial and final configurations are the same. 
It is important to remember that drawing paths in the plane is somewhat misleading: as 
in three dimensions, opposite points are identified, so that a path from any point to the 
diametrically opposite point is closed. 




Figure 8: Three topologically distinct paths in two dimensions. 



The first and second paths are similar to the direct and exchange paths of the three- 
dimensional problem. The third path, however, represents a distinct class of path in two 
dimensions. The particles circle around each other, returning to their starting points. 

It is perhaps easier to visualize these paths in a three-dimensional space-time plot, where 
the vertical axis represents time and the horizontal axes represent space (Figure 9). 

It is clear that in the third path the particle initially at qi returns to qi, and similarly 
for the other particle: this path does not involve a permutation of the particles. It is also 
clear that this path cannot be continuously deformed into the first path, so it is in a distinct 
topological class. (It is critical here that we have excised the origin in relative coordinates - 
i.e., that we have disallowed configurations where the two particles are at the same point in 
space.) 
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The existence of this third class of paths generahzes in an obvious way, and we are led 
to the following conclusion: the paths starting and finishing at relative position q can be 
divided into an infinite number of classes of paths in the plane (minus the origin); a class is 
specified by the number of interchanges of the particles (keeping track of the sense of each 
interchange). This is profoundly different from the three-dimensional case, where there were 
only two classes of paths: direct and exchange. 

If we characterize a path by the polar angle of the relative coordinate, this angle is mr 
in the n*^ class, where n is an integer. (For the three paths shown above, n = 0, 1, and 2, 
respectively.) 

We can write 

oo 

A(q,T;q,0)= ^ ^(q,T;q,0), 

n=— oo 

where An is the covering- space PI considering only paths of change of polar angle mr. 
This path integral can again be generalized to 

oo 

^(q,r;q,0)= ^ C„yl„(q, T; q, 0), 

n=— oo 

Cn being phases. Since each An satisfies the Schroedinger equation, so does this generaliza- 
tion. 

Again, a restriction on the phases arises, as can be seen by the following argument. Let 
us relax the condition that the initial and final points are the same; let us denote the final 
point q' by its polar coordinates (g', 6'). Writing A{q', 9') = A{c^, T; q, 0), we have 

oo 

A{q',9') = Y.CnAn{q',9'). 

— oo 

Now, we can change continuously 9' —>■ 9' + i:, yielding 

oo 

A{q',9' + 7r)= J2 CnAniq',9' + 7r). (26) 

n=— oo 




Figure 9: Spacetime depiction of the 3 paths in Figure 8. 
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Two critical observations can now be made. First, the final configuration is unchanged, so 
A[q', 6' + tt) can differ from A{q', 9') by at most a phase: 

A(g',^' + 7r) = e-*^A(g',^')- 

Second, A„(g', 6' + ti) = An+i{q', 9'), since this is just two different ways of expressing exactly 
the same quantity. Applying these two observations to (|26|) , 

oo oo 

e-*-^ ^ CnAniq',9') = Y.CnAn+iiq',9') 

n=—oo — oo 

oo 

= Y.Cn-iAn{q',9'). (27) 

— oo 

Equating coefficients of An{q',9')^ we get 
Choosing Co = 1, we obtain for the amplitude 

oo 

A = ^e^"^^, (28) 

— oo 

which is the two-dimensional analog of the three-dimensional result (PTf). 

The most important observation to be made is that there is no longer a restriction on 
the angle 0, as was the case in three dimensions. We see that, relative to the "naive" PI 
(that with = 0), the class corresponding to a net number n of counter-clockwise rotations 
of one particle around the other contributes with an extra phase expin(f). If = or tt, 
this collapses to the usual cases of Bose and Fermi statistics, respectively. However in the 
general case the phase relation between different paths is more complicated (not determined 
by whether the path is "direct" or "exchange"); this new possibility is known as fractional 
statistics, and particles obeying these statistics are known as anyons. 

Anyons figure prominently in the accepted theory of the fractional quantum Hall effect, 
and were proposed as being relevant to high-temperature superconductivity, although that 
possibility seems not to be borne out by experimental results. Perhaps Nature has other 
applications of fractional statistics which await discovery. 

4.3 Magnetic Monopoles and Charge Quantization 

All experimental evidence so far tells us that all particles have electric charges which are 
integer multiples of a fundamental unit of electric charge, e.Q There is absolutely nothing 
wrong with a theory of electrodynamics of particles of arbitrary charges: we could have 
particles of charge e and a/TTc, for example. It was a great mystery why charge was quantized 
in the early days of quantum mechanics. 

In 1931, Dirac showed that the quantum mechanics of charged particles in the presence of 
magnetic monopoles is problematic, unless the product of the electric and magnetic charges 

^The unit of electric charge is more properly e/3, that of the quarks; for simplicity, I will ignore this fact. 
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is an integer multiple of a given fundamental value. Thus, the existence of monopoles 

implies quantization of electric charge, a fact which has fueled experimental searches for and 
theoretical speculations about magnetic monopoles ever since. 

We will now recast Dirac's argument into a modern form in terms of Pis. A monopole 
of charge g positioned at the origin has magnetic field 




As is well known, this field cannot be described by a normal (smooth, single-valued) vector 
potential: writing B = V x A implies that the magnetic flux emerging from any closed 
surface (and thus the magnetic charge contained in any such surface) must be zero. This 
fact makes hfe difficult for monopole physics, for several reasons. Although classically the 
Maxwell equations and the Lorentz force equation form a complete set of equations for 
particles (both electrically and magnetically charged, with the simple addition of magnetic 
source terms) and electromagnetic fields, their derivation from an action principle requires 
that the electromagnetic field be described in terms of the electromagnetic potential, A^. 
Quantum mechanically, things are even more severe: one cannot avoid ^4^, because the 
coupling of a particle to the electromagnetic field is written in terms of A^, not the electric 
and magnetic fields. 

Dirac suggested that if a monopole exists, it could be described by an infinitely-thin and 
tightly-wound solenoid carrying a magnetic flux equal to that of the monopole. The solenoid 
is semi-infinite in length, running from the position of the monopole to infinity along an 
arbitrary path. The magnetic field produced by such a solenoid can be shown to be that of 
a monopole plus the usual field produced by a solenoid, in this case, an infinitely intense, 
infinitely narrow tube of flux running from the monopole to inflnity along the position of 
the solenoid (Figure 10). Thus except inside the solenoid, the held produced is that of the 
monopole. The held inside the solenoid is known as the "Dirac string" . The flux brought 




Figure 10: Monopole as represented by a semi-infinite, infinitely tightly-wound solenoid. 
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into any closed surface including the monopole is now zero, because the solenoid brings in a 
flux equal to that flowing out due to the monopole. Thus, the combined monopole-solenoid 
can be described by a vector potential. 

However, in order for this to be a valid description of the monopole, we must somehow 
convince ourselves that the solenoid can be made invisible to any electrically charged particle 
passing by it. We can describe the motion of such a particle by a PI, and two paths passing 
on either side of the Dirac string can contribute to the PI (Figure 11). But the vector 
potential of the Dirac string will affect the action of each of these paths differently, as we 
have seen in the Aharonov-Bohm effect. 




Figure 11: Paths contributing to the propagator in the presence of a monopole. The paths 
form a loop encircling the Dirac string. 

In order that the interference of these paths be unaffected by the presence of the Dirac 
string, the relative phase must be an integral multiple of 27r. This phase is 

e / e$ ineg 



(j) dq - 



he J he he 

Setting this to 27rn, in order for the motion of a particle of charge e to be unaffected by the 
presence of the Dirac string, the electric charge must be 

2nhe he 

e = n = —n. (29) 

Thus, the existence of magnetic monopoles requires the quantization of electric charge; the 
fundamental unit of electric charge is 27ihe/g. 

In modern theories of fundamental physics. Grand Unified Theories also imply quanti- 
zation of electric charge, apparently avoiding the necessity for magnetic monopoles. But 
any Grand Unified Theory actually has magnetic monopoles as well (though they are of a 
nature quite different to the "Dirac monopole"), so the intimate relation between magnetic 
monopoles and the quantization of electric charge is preserved, albeit in a form quite different 
from that suggested by Dirac. 



5 Statistical Mechanics via Path Integrals 

The path integral turns out to provide an elegant way of doing statistical mechanics. The 
reason for this is that, as we will see, the central object in statistical mechanics, the partition 
function, can be written as a PI. Many books have been written on statistical mechanics 
with emphasis on path integrals, and the objective in this lecture is simply to see the relation 
between the partition function and the PI. 
The definition of the partition function is 

Z = ^e-^^^ (30) 
j 

where (3 = 1/kBT and Ej is the energy of the state We can write 

j 

But recall the definition of the propagator: 

W,T;g,0) = (g'|e-'^^|g). 

Suppose we consider T to be a complex parameter, and consider it to be pure imaginary, so 
that we can write T = where (3 is real. Then 

= (9le-^^Eb)0ii?) 

j 

=1 

Putting q' = q and integrating over q, we get 

ldqK{q,-tf3;q,0)=J2e-^'^^{j\ J dq\q) {q\ \j) = Z. (31) 

This is the central observation of this section: that the propagator evaluated at negative 
imaginary time is related to the partition function. 

We can easily work out an elementary example such as the harmonic oscillator. Recall 
the path integral for it, ([T7|): 

We can put q' = q and T = — i/5: 

/ mtu \^^^ f mujq'^ , /n \ A 
K{q, -1(3; g, = - — . , , exp <^ — . x cosh " ^ f • 



5 STATISTICAL MECHANICS 



26 



The partition function is thus 



Z ^ J dqK{q,-if3;q,0)^(^ 



1/2 



TT 



[2(cosh(/3cu) - 1)] 



-1/2 



27r sinh(/9c<;) y 



sinh(/3aj) 



(cosh (/3a;) 



-Pui/2 



-/30-+l/2)a; 



1 - e->^^ 

Putting h back in, we get the famihar result 

CXD 

Z = j:e 

j=0 



-/3(i+l/2)na; 



The previous calculation actually had nothing to do with Pis. The result for K was 
derived via Pis earlier, but it can be derived (more easily, in fact) in ordinary quantum 
mechaincs. However we can rewrite the partition function in terms of a PI. In ordinary 
(real) time. 



K{q', T- q:0)^J Vq{t) exp ^ ^ dti^- V{q) j , 



where the integral is over all paths from (g, 0) to (g', T). With q' = q, T 



K{q, -i/3;g,0)= Jvq{t)expiJ^ 



dt — f- 



V{q) . 



where we now integrate along the negative imaginary time axis (Figure 12). 
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Figure 12: Path in the complex time plane. 

Let us define a real variable for this integration, t — it. r is called the imaginary time, 
since when the time t is imaginary, r is real. (Kind of confusing, admittedly, but true.) Then 

the integral over r is along its real axis: when t : ^ —iP, then r : ^ j3. Wc can write q 
as a function of the variable r: q{t) —>■ q{T); then q = idq/dr. The propagator becomes 



K{q, -1(3- q,0) = J Vq{T) exp - £ dr ( ) + V{q) 



(32) 
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The integral is over all functions g(r) such that g(0) = = q. 

The result ( ^2|) is an "imaginary-time" or "Euclidean" path integral, defined by asso- 
ciating to each path an amplitude (statistical weight) exp —Se, where Se is the so-called 
Euclidean action, obtained from the usual ( "Minkowski" ) action by changing the sign of the 
potential energy term. 

The Euclidean PI might seem like a strange, unphysical beast, but it actually has many 
uses. One will be discussed in the next section, where use will be made of the fact that at low 
temperatures the ground state gives the dominant contribution to the partition function. It 
can therefore be used to find the ground state energy. We will also see the Euclidean PI in 
Section 9, when discussing the subject of instantons, which are used to describe phenomena 
such as quantum mechanical tunneling. 



6 Perturbation Theory in Quantum Mechanics 



We can use the Euclidean PI to compute a perturbation expansion for the ground state energy 
(among other things). This is not terribly useful in and of itself (once again, conventional 
techniques are a good deal easier) , but the techniques used are very similar to those used in 
perturbation theory and Fcynman diagrams in field theory. For this reason, we will discuss 
corrections to the ground state energy of an elementary quantum mechanical system in some 
detail. 

Prom Z it is quite easy to extract the ground state energy. (This is a well-known fact of 
statistical mechanics, quite independent of Pis.) Prom the definition of Z, 



we can see that the contribution of each state decreases exponentially with /?. However, that 
of the ground state decreases less slowly than any other state. So in the hmit of large /3 {i.e., 
low temperature), the ground state contribution will dominate. (This is mathematically 
straightforward, and also physically reasonable.) One finds 

£;o = - lim ^logZ. (33) 

In fact, we can extract Eq from something slightly easier to calculate than Z. Rather 
than integrating over the initial (= final) position, as with Z, let us look at the Euclidean 
propagator from g = to g' = (the choice of zero is arbitrary). 

XE(0,/3;0,0) = (g' = 0|e-^^|g = 0). 

We can insert a complete set of eigenstates of H: 

X^(0,/3;0,0) = (g' = 0|e-^^$:ii)(i|g = 0) 

3 

where 4>j are the wave functions of H. Again the ground state dominates as /? — > oo, and 

£;o = - lim ^logXs(0,/3;0,0). (34) 

(As j3 — >■ cxo, the difference between f3~^\ogZ and (3~'^ log Ke goes to zero.) 

So let us see how we can calculate Ke{0, jS; 0, 0) perturbatively via the PI. The starting 
point is 

where the paths over which we integrate start and finish at g = 0, and where the Euclidean 
action is 
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and with q = dq/dr. As an example, consider the anharmonic oscillator, with quadratic and 
quartic terms in the potential: 



Ke{0, f3; 0, 0) = 1 Vq exp - J dr (^^mq^ + ^muj\^ + Ag^^ . 



(35) 



Clearly it is the quartic term which complicates life considerably; we cannot do the PI 
exactly.!] But we can use the following trick to evaluate it perturbatively in A. (This trick is 
far more complicated than necessary for this problem, but is a standard - and necessary! - 
trick in quantum field theory.) Define i^£;[J], the PI for a harmonic oscillator with a source 
term (which describes the action of an external force) added to the Lagrangian: 

[J] = JVqexp- J dr Qmg^ + ^mu^q"^ - J(r)g(r)) . (36) 

Unlike (|35D, this PI can be evaluated exactly; we will do this (as much as is necessary, at 
least) shortly. Once we have evaluated it, how does it help us to compute (^5])? To see the 
use of K^[J], acting on it with a derivative has the effect of putting a factor q in the PI: for 
any time Ti, 

-j-jj—^ = JVqq{n)exp-JdT\^-mq + -muj q -J{T)q{T) 
A second derivative puts a second q in the PI: 

6Jin)SJir,) = J ^^g("^)^("^)^"P-y dr [-mq + -rm. q - J{r)q{ 
In fact, we can generalize this to an arbitrary functional F: 

Kl[J]= jVqF[q]e-''-^-^\ (37) 



5J 



where S^yJ] is the Euclidean action for the harmonic oscillator with source. (To prove (P7|), 
bring F[5/5J] inside the PI; each 5/5J in F operating on exp— ^^[J] gives rise to a g in 
front of the exponential.) 

Now, if we choose F[q\ = exp — / dr^q^, we get: 

e-/'^-^(^)V°[J] = |Pgexp|-|rfr^g4|e-50[J] 

= Jvqexp- J dr (^mq^ + ^mu^q^ + - J (r)g(r) j . 



^In fact, the situation is exactly like the evaluation of the ordinary integral 

/°° 1 A 

4 4! 

which looks innocent enough but which cannot be evaluated exactly. The technique which we will develop 
to evaluate (pq) can also be used for this ordinary integral - an amusing and recommended exercise. 
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If we now put J = 0, we have the PI we started with. So the final result is: 



KEiO,(3;0,0)= [expi- J 



J=0 



(38) 



We can, and will, calculate as an explicit functional of J. If we then expand the 

exponential which operates on it in (pSj), we get a power series in A: 



dr 



X 



4! \SJ{t) 



6 



dr' 



X 



6 



4! V^^(^'), 



+ 



j=o 



oiX' 



J=0 



Let us now evaluate K'^lJ], 

K'^[J] = J Vq exp — J dr (^-mq'^ + -muj'^q'^ — J{T)q{T) 

To do this, suppose that we can find the classical path qcj{T), the solution of 

mq = muj'^q — J{t), g(0) = g(/?) = 0. 



(39) 



Once we have done this, we can perform a change of variables in the PI: we define g(r) 
9c J (t) + 2/(t), and integrate over paths yir). This is useful because 



j dr Qmg^ + ^mu'^q^ - J{T)q{T)j = J dr (^mqlj + ]^muj'^qlj - J{T)q^j{T)j 

+ (linear in y) + j dr (^-my'^ + -mu'^y'^^ . 

=0 

The linear term vanishes because qcj satisfies the equation of motion. So the PI becomes 



KM 



-Sec[J] 



Vy exp — j dr ^-my'^ + -mu'^y'^^ . 



The crucial observation is that the resulting PI is independent of J: it is an irrelevant 
constant; call it C. (In fact, C is neither constant [it depends on nor entirely irrelevant 
[it is related to the unperturbed ground state energy, as we will see] . Crucial for the present 
purposes is that C is independent of J.) 



where 



Sec[J] = j dr (^^mqlj + ^mu^qlj - J{T)q^j{r] 
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This can be simplified by integrating the first term by parts, yielding 

Sec[J] = J drq^j (^-^mqcj + ^mu^qcj - J(- 



Using the classical equation of motion (|39|), we get 



We must still solve the classical problem, (p9|). The solution can be written in terms of 
the Green's function for the problem. Let G{t, t') be the solution of 



m 



-c.MG(r,r') = 5(r-r'), 



G(0,r') = G(/3,r') 



Then we can immediately write 

qcj{r)= f dr'G{r,T')J{T'), 
Jo 

which can be proven by substution into (pQl). We can now write 



Kl[J] = C exp^J drdr' J{t)G{t, t')J{t'). 



(40) 



We can find the Green's function easily in the limit (3 ^ oo. It is slightly more convenient 
to treat the initial and final times more symmetrically, so let us choose the time interval to 
be (— /3/2, +P/2); in the limit /? ^ oo we go from — oo to oo. Then we have 

m(^-u;A G{T,T') = 6{T-r'). 



By taking the Fourier transform, we see that 



^ik(T-T') 



m J-oo 271 [hp' + uj"^ 



(41) 



We can now compute the first-order correction to Ke (from which we get the first-order 
correction to the ground state energy). We have 



(42) 



J=0 



Since in the second term we take four derivatives of and then set J = 0, only the 

piece of which is quartic in J is relevant: fewer than four J's will be killed by the 
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derivatives; more than four will be killed when setting J = 0. 

K° [J] = C exp^J drdr' J{t)G{t, t')J{t') 

= irrelevant + ^ " ^ / drdr' J{t)G{t, r') J(r )) 

= ^{JiGuJ2){J3G3,J4), (43) 



where we have used the compact notation {J1G12J2) = J dTidT2 J{ti)G{ti,T2)J{t2) 
Substituting (P|) into (P), 



KE = c(^l-^ydT (^^j {JiGuJ2){J3GuJ4) + o{X'yj . (44) 

To ensure that we understand the notation and how functional differentiation works, let 
us work out a slightly simpler example than the above. Consider 



X = 



f \ 2 / r \ 2 





6Jir) 



J1G12J2) = [jj^^ J dndT2Jin)Gin,r2)J{T2). 



The first derivative can act either on Ji or J2. In either case, it gives a delta function, which 
will make one of the integrals collapse: 



X = J dndT2 (5(r - n)G{n, ra) J(r2) + J{n)G{n, r)5(r - rs)) 

^ ^ J dT2G{T,T2)J{T2) + J rfn J(ri)G(ri , r) 



6J{r) 

In each term the remaining derivative acts similarly and kills the remaining integral; the 
result is 

X = 2G{t,t). 

The functional derivatives in (^31) are a straightforward generalization of this; we find 

KE = G(^l-iydT4\G{r,rfy 
From (|41|) G{t, t) = —l/2mu] the r integral is just the time interval /?, and finally we get 
Ke{0, P; 0, 0) = C (^1 - + o{\')^ = Ce-^V32n^^.^ ^45^ 

to order A. 

Now we can put this expression to good use, extracting the ground state energy from 

my- 

Eo = - lim ^ log KEiO, /3; 0,0) = - lim ^ (log C * 
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Recall that the constant C in (|45|) depends on (3; this dependence must account for the 
ground state energy; the term linear in A gives the first correction to the energy. Thus 

where we have reintroduced h. We can check this result against standard perturbation theory 
(which is considerably easier!); the first-order correction to the ground state energy is 



as above. 

One's sanity would be called into question were it suggested that the PI calculation is a 
serious competitor for standard perturbation theory, although the latter itself gets rapidly 
more and more messy at higher orders. The technique above also gets messier, but it may 
well be that its messiness increases less quickly than that of standard perturbation theory. If 
so, the PI calculation could become competitive with standard perturbation theory at higher 
orders. But really the main motivation for discussing the above method is that it mimics in 
a more familiar setting standard perturbation techniques in quantum field theory. 

To summarize this long and somewhat technical section, let us recall the main features 
of the above method. We express the ground state energy as an expression involving the 
propagator, (pi]). We separate the Lagrangian into a "free" {i.e., quadratic) part and an 
"interacting" (beyond quadratic) part. Via the PI, we write the interacting propagator in 
terms of a free propagator with source term added, (pHj); this expression is amenable to a 
perturbation expansion. The free propagator can be evaluated explicitly, (^OD ; then ( ^8]) can 
be computed to any desired order. From this we obtain directly the ground state energy to 
the same order. 



7 Green's Functions in Quantum Mechanics 



T(^(xi)0(x2) 



In quantum field theory we are interested in objects such as 

(O|r0(Xi)0(x2)---0(x„)|O), 

the vacuum expectation value of a time-ordered product of Heisenberg field operators. This 
object is known as a Green's function or as a correlation function. The order of the operators 
is such that the earliest field is written last (right-most), the second earliest second last, etc. 
For example, 

J 0(Xi)0(x2) > X2 

\ 0(x2)0(Xi) X^2 > ^1 

Green's functions are related to amplitudes for physical processes such as scattering and 
decay processes. (This point is explained in most quantum field theory books.) 
Let us look at the analogous object in quantum mechanics: 

G^^\h, t2, ■ ■ ■ , tn) = (0| Tq{t,)q{t2) ■ ■ ■ q{tn) |0) . 

We will develop a PI expression for this. 

First, we must recast the PI in terms of Heisenberg representation objects. The operator 
q{t) is the usual Heisenberg operator, defined in terms of the Schroedinger operator q by 
q{t) = e^^^qe^'^^K The eigenstates of the Heisenberg operator are \q,t): q(t) \q,t) = q \q,t). 
The relation with the time-independent eigenstates is \q,t) = e*^* Then we can write 
the PI: 

K={q'\ e-'^^ \q) = {q', T\ g, 0) = J Vq e'^ . 

We can now calculate the "2-point function" (^(^1,^2), via the PI. We will proceed in two 
steps. First, we will calculate the following expression: 

(g',T|Tg(ti)g(t2)|g,0). 

We will then devise a method for extracting the vacuum contribution to the initial and final 
states. 

Suppose first that ti > ^2- Then 

{q\T\Tq{h)q{h)\q,Q) = {q' ,T\ qitMh) hO) 

= / dqidq2{q',T\qi,ti) {qi,ti\q{ti)q{t2)\q2,t2) {q2,t2\q,0) 
J V V " 

(gi,ti|gi q2\g2,t2) 



J dqidq2qiq2{q,T\qi,ti) (gi,ti| 52,^2) (^2,^2! g, 0). 



^There is a possible point of confusion here. We all know that Heisenberg states are independent of time, 
yet the eigenstates of q{t) depend on time. Perhaps the best way to view these states \q,t) is that they 
form, for any fixed time, a complete set of states. Just like the usual (time-independent) Heisenberg state 
\q) describes a particle which is localized at the point q at time t — 0, the state \q,t) describes a particle 
which is localized at the point q at time t. 
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Each of these matrix elements is a PI: 

dqidq2qiq2 / Vqe''' / Vqe'^ / Pge*^. 

This expression consists of a first PI from the initial position q to an arbitrary position q2, 
a second one from there to a second arbitrary position qi , and a third one from there to the 
final position q'. So we are integrating over all paths from q to q', subject to the restriction 
that the paths pass through the intermediate points qi and q2- We then integrate over the 
two arbitrary positions, so that in fact we arc integrating over all paths: we can combine 
these three path integrals plus the integrations over qi and q2 into one PI. The factors qi 
and q2 in the above integral can be incorporated into this PI by simply including a factor 
q{ti)q{t2) in the PI. So 

{q',T\q{h)q{t2) \q,0) = Vqq{tMh)e'' {h > ^2). 

Jg,0 

An identical calculation shows that exactly this same final expression is also valid for 
t2 < h: magically, the PI does the time ordering automatically. Thus for all times 

{q',T\Tq{h)q{t2)\q,0) = ["'^ Vqq{h)q{t2)e'^ . 

Jq,0 

As for how to obtain vacuum-to- vacuum matrix elements, our work on statistical mechan- 
ics provides us with a clue. We can expand the states (g', T| and \q, 0) in terms of eigenstates 
of the Hamiltonian. If we evolve towards a negative imaginary time, the contribution of all 
other states will decay away relative to that of the ground state. We have (resetting the 
initial time to — T for convenience) 

{q',T\q, -r)oc(o,r|o,-r), 

where on the right the "0" denotes the ground state. The proportionahty involves the ground 
state wave function and an exponential factor exp2iEoT = exp — 2£'o|T|. 

Wc could perform all calculations in a Euclidean theory and analytically continue to real 
time when computing physical quantities (many books do this), but to be closer to physics 
we can also consider T not to be pure imaginary and negative, but to have a small negative 
imaginary phase: T — |T|e~*^ (e > 0). In what follows, I will simply write T , hut please keep 
in mind that it has a negative imaginary part! With this, 

(0,T|0,-T)oc(g',T|g, -T) = j Vqe'^ . 

To compute the Green's functions, we must simply add Tq{ti)q{t2) ■ ■ -qitn) to the matrix 
element, and the corresponding factor q{t\)q{t2) • ■ - qitn) inside the PI: 

(0, T\ Tq{t^)q{t2) ■ ■ ■ q{tn) |0, -T) ^ j Vq q{ti)q{t2) ■ ■ ■ q{tn)e'^ ■ 
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The proportionality sign is a bit awkward; fortunately, we can rid ourselves of it. To do this, 

wc note that the left hand expression is not exactly what we want: the vacua |0, ±T) differ 
by a phase. We wish to eliminate this phase; to this end, the Green's functions are defined 

G^^\h,t2,---,tn) = {0\Tq{t,)q{h)---q{tn)\0) 

^ {0,T\Tq{ti)qih)---qitn)\0,-T) 

{0,T\0,-T) 

!Vqq{h)q{t2)---q{tn)e'^ 
JVqe^s 

with no proportionality sign. The wave functions and exponential factors in the numerator 

and denominator cancel. 

To compute the numerator, we can once again use the trick we used in perturbation 
theory in quantum mechanics, namely, adding a source to the action. We define 



Z[J] 



JVqe^s 

If we operate on Z[J] with i~^S/5J{ti), this gives 



(0|0) 



1 S 



i 5J{ti) 



Z[J] 



JVqe^^ 

fVqq{h)e'^ 

JVqe^ 
(0,T|g(ti)|0,-T) 



j=o 



(0,T|0,-T) 



= (0| m |0) 



(The expectation values are evahiated in the absence of J.) 

Repeating this procedure, we obtain a PI with several g's in the numerator. This ordinary 
product of g's in the PI corresponds, as discussed earlier in this section, to a time-ordered 
product in the matrix element. So we make the following conclusion: 



i 5J{ti) i 5J{tr. 



-Z[J] 



jVqq{ti) ■ ■■q{tn)e 



iS 



J=0 



jVqe 



iS 



Tq{t,)---q{t,)\0). 



For obvious reasons, the functional Z[J] is called the generating functional ior Green's func- 
tions; it is a very handy tool in quantum field theory and in statistical mechanics. 
How do we calculate Z[J]7 Let us examine the numerator: 

N = jVqe'^^^-f 'itJ{t)q{t))_ 
Suppose initially that -S" is the harmonic oscillator action (denoted -S'o): 

Sq^ J dt Qmg^ - ^muj'^q'^^ , 
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Then the corresponding numerator, Nq, is the non-Euchdean [i.e., real-time) version of the 
propagator -^'^'[J] we used in Section 6. We can calculate No[J] in the same way as K'^lJ]. 
Since the calculation repeats much of that of /r£;[J], we will be succinct. 
By definition, 



No = J Vq{t) expz J dt Qmg^ 



We do the path integral over a new variable y, defined by qit) = qdt) + y{t), where qc is 
the classical solution. Then the PI over ?/ is a constant (independent of J) and we can avoid 
calculating it. (It will cancel against the denominator in Z[J].) Calling it C, we have 



where 



SQ.j[qc 



dt 



]^mql - ^muj'^ql + Jq^ 



dtJ{t)q,{t), 



using the fact that qc satisfies the equation of motion. We can write the classical path in 
terms of the Green's function (to be determined shortly), defined by 



G{t,t') = -i5{t-t'). 



(46) 



Then 



We can now write 



q^it) = -i / dt'G{t,t')J{t'). 



No = Cexp^J dtdt' J{t)G{t, t')J{t') 



Dividing by the denominator merely cancels the factor C, giving our final result: 



Z[J] = expiy" dtdt' J{t)G{t,t')J{t'). 



We can solve (|46| ) for the Green's function by going into momentum space; the result is 

dk i 



G{t,t') =G{t-t') 



-ik(t-t') 



27r P - 

However, there are poles on the axis of integration. (This problem did not arise in Euclidean 
space; see (^Tf) .) The Green's function is ambiguous until we give it a "pole prescription", 
i.e., a boundary condition. But remember that our time T has a small, negative imaginary 
part. We require that G go to zero as T — oo. The correct pole prescription then turns out 
to be 



G{t - 1') 



dk 



2tt k'^ — uj'^ + ie 



-ik(t-t') 



(47) 



We could at this point do a couple of practice calculations to get used to this formalism. 
Examples would be to compute perturbatively the generating functional for an action which 
has terms beyond quadratic (for example, a term), or to compute some Green's function in 
either the quadratic or quartic theory. But since these objects aren't really useful in quantum 
mechanics, without further delay we will go directly to the case of interest: quantum field 
theory. 



8 Green's Functions in Quantum Field Theory 



It is easy to generalize the PI to many degrees of freedom; we have in fact already done so in 
Section 4, where particles move in two or three dimensions. It is simply a matter of adding 
a new index to denote the different degrees of freedom (be they the different coordinates 
of a single particle in more than one dimension or the particle index for a system of many 
particles). 

One of the most important examples of a system with many degrees of freedom is a field 
theory: q{t) — > 0(x, t) = Not only is this a system of many degrees of freedom, but 

one of a continuum of degrees of freedom. The passage from a discrete to continuous system 
in path integrals can be done in the same way as in ordinary classical field theory: we can 
discretize the field (modeling it by a set of masses and springs, for instance), do the usual 
path integral manipulations on the discrete system, and take the continuum limit at the 
end of the calculation. The final result is a fairly obvious generalization of the one-particle 
results, so I will not dwell on the mundane details of discretization and subsequent taking 
of the continuum limit. 

The analog of the quantum mechanical propagator is the transition amplitude to go from 
one field configuration 0(x) at t = to another 0'(x') at t = T: 

i^(0'(x'),T;0(x),O) = J P0e*^[*l, (48) 
where S is the field action, for instance 

S[<P]= / rf^xQ (9,0)2 -ImV) (49) 

for the free scalar field. In (^) the integral is over all field configurations 0(x) obeying the 
stated initial and final conditions. 

In field theory, we are not really interested in this object. Rather (as mentioned earlier), 
we are interested in Green's functions. Most of the work required to translate (p8[) into an 
expression for a Green's function (generating functional of Green's functions, more precisely) 
has already been done in the last section, so let us study a couple of cases. 

8.1 Free scalar field. 

For the free scalar field, whose action is given by (|49|) , the generating functional is 

(0|0)j=„ 

Both numerator and denominator can be written in terms of Pis. The numerator is 
No = jVct> expz j d^x {^{d^ct>f - ^m^^^ + j^^ 

We write = 0c + V^, where 0c is the classical field configuration, and integrate over the 
deviation from 0c. The action can be written 



^[0c + v] = j d'x (^(9,0c)2 - ^m'<Pl + J0c) + / d'x (^^{d.^r - \ 



2 2 

m if 



5 
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where as usual there is no term hnear in (p since (f)c by definition extremizes the classical 
action. So 

No = Cexpt J d'x - lm^(pl + J(p?j, 



where 



C is independent of J and will cancel in Z. (Indeed, the denominator is equal to C] 
Using the fact that 0c obeys the classical equation 

{d^ + m2)0, = J, 



we can write 



Nq = C exp - / d'^x J{x)(f)c{x). 



Finally, we can write 0c in terms of the Klein-Gordon Green's function, defined by 

{d^ + m^)AF{x, x') = -i5\x - x'). 



It is 



so 



(f)cix) = i d'^x Ap{x, x')J{x'), 



Zq = — ^ = exp — / d'^xd'^x' J{x)Ap{x,x')J{x') 
C 2 J 



The Green's function is found by solving its equation in 4-momentum space; the result is 

d^k i 



Af{x,x') 



(27r)4 P - m2 + ie 



e-*^-(^-^') = Ap(x - 



X 



adopting the same pole prescription as in p7|). Note that Ai? is an even function, Af{x — 
x') = Af{x' — x). 

Let us calculate a couple of Green's functions. These calculations are reminiscent of those 
at the end of Section 6. As a first example, consider 



Gj^\xi,X2) = {0\T^{xi)4>{x2) |0) 



-Zo[J] 



J=0 



P \6J{xi)6J{x2) 
Expanding Zq in powers of J, 

Zo[J] = l-^J dVx' J{x)Af{x - x')J{x') + o( J^). 

The term quadratic in J is the only one that survives both differentiation (which kills the 
"1") and the setting of J to zero (which kills all higher-order terms). So 



Gq\xi,X2) 



P5J{xi)5J{x2) 



J d^xd^x' J{x)Af{x - x')J{x') 
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There arise two identical terms, depending on which derivative acts on which J. The result 
is 

(2) 

Go '{xi,X2) = Af{xi - X2). 

So the Green's function (or two-point function) in the quantum field theory sense is also the 
Green's function in the usual differential-equations sense. 
As a second example, the four-point Green's function is 

G'o\xi,X2,X3,X4) = ^ ( — ^— • • ■ — |— exp-i / d'^xd'^x' J{x)Af{x,x')J{x') 

\oJ[x\) dJ{x4) 2 J 

This time the only part of the exponential that contributes is the term with four J's. 



j=o 



G^^\X,,X2,X„X,) 



d xd x' J{x)Af{x,x')J{x'] 



6J{xi) SJ{x4) 2 V 2 

There are 4! = 24 terms, corresponding to the number of ways of associating the derivatives 
with the J's. In 8 of them, the Green's functions which arise are Af{xi — X2)Af{x3 — X4), 
and so on. The result is 

G''q\xi,X2,X3,X4) = Af{xi - X2)Af{x3 - X4) + Af{xi - X3)Af{x2 - X4) 

+Af{xi - X4)Af{x2 - X3), (50) 

which can be represented diagramatically as in Figure 13. 



G 



(4) 




Figure 13: Diagrammatic representation of (^Ol). Each line counts as a factor with 
argument corresponding to the endpoints of the line. 



8.2 Interacting scalar field theory. 

Usually, if the Lagrangian has a term beyond quadratic we can no longer evaluate exactly the 
functional integral, and we must resort to perturbation theory. The generating functional 
method is tailor-made to do this in a systematic fashion. To be specific, consider 0^ theory, 
defined by the Lagrangian density 

Then the generating functional is (up to an unimportant constant: we will normalize ulti- 
mately so that Z[J = 0] = 1) 

Z[J] =cjv<pexptj d'x (l{d,<P)' - ^m'<p' - ^<p' + J A 
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Because of the quartic term, we cannot evaluate the functional integral exactly. But we 

can use a trick wc first saw when discussing perturbation theory in quantum mechanics: 
replacing the higher-order term by a functional derivative with respect to J: 

Z[J] = C jV(j) exp j expi j d^x - ^m'^^ + J^) 

We can pull the first exponential out of the integral; the functional integral which remains 
is that for Zq. Adjusting the constant C so that Z[J = 0] = 0, we get 

exp j— i^/ d'^x (^jjp-^ Icxp— |/ d'^xd^x' J{x)Ap{x, x')J{x') 

z[j] = ^ ■ ^ 



1 V(j) exp < 




[-4.1 






y Vcj) exp < 


[-4J 



^exp|-iA/ci%(i^) |exp-i/ci%(iVJ(x)A^(x,x')J(x')^ 

This expression now enables us to compute a perturbative expansion for any Green's function 
we desire. This is a rather mechanical job, and the only way to learn it is by doing lots of 
examples. To illustrate the method, let us look at G^'^^xi, X2) to the first nontrivial order in 
A. 

We have 



G^'\X„X2) 



()./(,) 


52 


) exp j 




1 () 


/(•'■)) J 


■cxp-i(./„A,v,,./,)^ 




J=0 




^exp < 


''■41 1 d X (j 




■ exp - 


-\{JaAFabJb)^ 


J=0 



where as in Section 6 (■ ■ ■) implies integration over the positions of the J's. In both numerator 
and denominator, wc can expand both exponentials. The only terms that survive are those 
that have the same total number of derivatives and J's. Let us look at the term linear in 
A in the numerator. There are six derivatives, so we need the term from the expansion of 
the second exponential with six J's. For this term, we get the following expression for the 
numerator: 

4 



5' 



5J{xi)5J{x2) 

There are now a total of 6! = 720 terms! However, only two distinct analytical expressions 
result. The first of these arises if the derivatives at Xi and X2 act on different (■ ■ ■)'s. A little 
combinatorial head-scratching tells us that there are 576 such terms, yielding the following 
expression: 

— —J d'^x Af{xi — x)Af{x — x)Af{x — X2), (51) 

which can be represented pictorially as in Figure 14 (a). 

The only other expression arises when the derivatives at xi and X2 act on the same (•••). 
This accounts for the remaining 144 terms; the analytic form which results is 

- ^-^Af{xi - X2) f d'^xApix - xf, (52) 
8 J 
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XXX X 

2 1 2 

(a) (b) 



Figure 14: Diagrammatic representation of (a) (|5TD and (b) (|5^) . Intersection points repre- 
sent vertices, and count as a factor —i\ J d^x. Symmetry factors complete the association of 
an analytic expression with each diagram. 



corresponding to the diagram in Figure 14(b). 

The denominator can be evaluated in a similar fashion; the Green's function to order A 



IS 



|Ai7(xi — X2) — y / d'^x Af{xi — x)Ap{x — x)Ap{x — X2) 

-f Ap(xi - X2) J d^x Apjx - xf + o(A^)} 

^""''""'^ 1-f /d%A^(x-a;)2 + o(A2) 

Since we have only computed the numerator and denominator to order A, we can rewrite 
this expression in the following way: 



G^'\X„X2) 



{Af{xi -X2) -Y-f d'^x^pixi - x)Af{x - x)Af{x - X2) + o(A^)} 
X {1 - fjd^xAFix - xf + o{X^)} 



1 - f J d^x Af{x - xy + o{X^) 



We can now cancel the second factor in the numerator against the denominator, to order A, 
resulting in 

z A f 

G'''^\xi, X2) = Af{xi — X2) — — j d'^x Af{xi — x)Af{x — x)Af{x — X2) + o(A^). 

This factorization of the numerator into a part containing no factors independent of the 
external position times the denominator occurs to all orders, as can be proven fairly cleanly 
via a combinatoric argument. The conclusion is that so-called disconnected parts (parts of 
diagrams not connected to any external line) cancel from Green's functions, a fact which 
simplifies greatly the calculation of these objects. 

It cannot be overemphasized that there are only three ways to get accustomed to this 
formalism: practice, practice, and practice. Other reasonable exercises are the calculation 
of G^^-* to order A^ and the calculation of G^'^^ to order A^. 0^ theory is also a useful testing 
ground for the techniques discussed in this section. 



9 Instantons in Quantum Mechanics 



9.1 General discussion 

It has already been briefly mentioned that in quantum mechanics certain aspects of a problem 
can be overlooked in a perturbative treatment. One example occurs if we have a harmonic 
oscillator with a cubic anharmonic term: V{q) = ^muj^q'^ + Xq^ (Figure 15). 




q 



Figure 15: V{q) = |ma;^g^ + Xq^. 

We can calculate corrections to harmonic oscillator wave functions and energies pertur- 
batively in A, to any desired order, blissfully ignorant of a serious pathology in the model. 
As can be seen from Figure 15, this model has no ground state: the potential energy is 
unbounded as g — — oo, a point completely invisible to perturbation theory. 

A second example is the double- well potential, V{q) = ^{q'^ — a^)^ (Figure 16). There 



V(q) 




-a a q 



Figure 16: V{q) = ^(g^ - a^- 

are two classical ground states. We can ignore this fact and expand V about one of the 
minima; it then takes the form of a harmonic oscillator about that minimum plus anharmonic 
terms (both cubic and quartic). We can then compute perturbative corrections to the wave 
functions and energies, and never see any evidence of the other minimum. Were we to expand 
about the other minimum, we would produce an identical set of perturbative corrections. By 
symmetry the ground state energies calculated perturbatively to any order will be the same 
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for the expansions about the two minima, so it appears that we have degenerate ground 
states. But in fact the ground state is not degenerate: a nonperturbative energy sphtting 
separates the true ground state (an even function of q) from the first excited state (an odd 
function); this sphtting is not seen in perturbation theory. 

We wiU examine this second example using Pis, the main goal being to calculate the 
energy splitting between the two candidate ground states. 

Let us first recall the PI expression for the Euclidean propagator: 



K 



2' 



-l3H/n 



Vq e 



-Ss/n 



where 



Se — 



f/3/2 / 1 \ 



Henceforth, we will set m 



1. Ke is useful because we can write it as 



(53) 



in the hmit /3 — > 00, this term will be dominated by the lowest-energy states. I say "states" 
here rather than "state" because we must calculate the two lowest-energy eigenvalues to 
get the splitting of the (perturbatively degenerate) lowest-energy states in the double-well 
potential. 

We will evaluate the PI using an approximation known as the semiclassical approxi- 
mation, or alternatively as the method of steepest descent. To illustrate it, consider the 
following integral 

J —00 

where S{x) is a function with several local minima (Figure 17). 



S(x) 




Figure 17: Potential with several minima. 

Suppose we are interested in this integral as — 0. Then the integral will be dominated 
by the minima of 5"; we can approximate it by a series of Gaussian integrations, one for each 
minimum of S. If Xi is such a minimum, then in its vicinity S{x) ~ S{xi) -|- |(x — XifS'^Xi)] 
we can write 

7~7i + 72 + 73 + ---, (54) 
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where 



/CO ^ 
dx exp —[S{xi) + -{x — XiYS"{xi)]/Ti 
-oo 2, 



Anharmonicities of S appear as corrections of order h to /. (This can be easily seen, for 
example, by considering a specific case such as S{x) = ax"^ + hx^.) 

We will compute the PI ([53|) in the semi-classical approximation, where the analog of the 
Xi in the above example will be classical paths (extremum of the action S'£;[g]). 

Suppose, then, that qdj) is the classical solution to the problem 

= qi-f3/2) = g, qiP/2) = q' . 

We can write g(r) = gc(T) + 1/(t); the action is 

r/3/2 



SE[qc + y] = (2^- + + ^(^- + ^) j 

= J ^dr {^-ql + V{qc)j + (linear in y) 



The term linear in y vanishes for the usual reason, and the higher order terms not written 
down are of cubic or higher order in y. Neglecting these (which give order h corrections to 
the PI), the propagator becomes 

Ke = J I)ge-^-/^ = e-^-[''d/"| Vyexp-J dr i^^f + ^V" {q,)y'^ /h. 

The functions ?/(r) over which we integrate satisfy the boundary conditions y{—(3/2) = 
y{(3/2) = 0. The PI, being Gaussian, can be done exactly; it is not as straightforward as the 
harmonic oscillator PI since V"{qc) depends on r. While we have often managed to avoid 
evaluating Pis, here we must evaluate it. (Unfortunately, this is rather difficult.) 

To this end, we can use a generalization of the Fourier expansion technique mentioned 
in Section 2.2.2. We can rewrite the action as 

SE = Jdr + ^l^"(ge)y^) =lldry + V"{q^)^ y. (56) 

The Schroedinger-like equation 

{--^ + V'iQc)) y = \y, yi-P/2) = yiP/2) = 
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has a complete, ortho normal set of solutions; let the solutions and eigenvalues be ?/fc(r) and 
Afc, respectively. The ort honor mality relation is 

13/2 

dTyk{T)yi{T) = 5ki. 

^13/2 

Then we can substitute ?/(r) = Y^k'^kUkiT) in (^6]), giving 

If ( I \ 1 [ 1 

Se = - J drY, akVk ( + V"{qc) 1 ^iVi = o ^ (^k^i^i J drykVi = o 51 ^l^k- 

k \ / I k,l k 

The PI can now be written as an integral over all possible values of the coefficients {ofc}. 
This gives 

Ke = J' !X{dake-^^<^''/^\ (57) 
k 

where J' is the Jacobian of the transformation from ?/(r) to {a^}. ( |57D is a product of 
uncoupled Gaussian integrals; the result is 

1 /2 

Ke = J'\[ (^) = J'm^^f'^W^k)-'" = J'n(2vr^)^/^det-V2 + V"{q^)\ , 

k \ / k k k \ / 

where we have written the product of eigenvalues as the determinant of the Schroedinger 
operator on the space of functions vanishing at ±/5/2. 
We can write J = J' nfc(27r^)^/^, giving 

Ke = Jdet-i/2 [-^ + ^"(^^)) (1 + ^(^))' 

where the o{h) corrections can in principle be computed from the neglected beyond-quadratic 
terms in (|53|). We will not be concerned with these corrections, and henceforth we will drop 
the {l + o{h)). 

9.2 Single Well in the Semiclassical Approximation 

Before looking at the double well, it is worthwhile examining the single well, defined be 

The classical equation is 

Note that this is the equation of motion for a particle moving in a potential —V{q). If we 
choose the initial and final points q = q' = 0, then the classical solution is simply qdT) = 0; 
furthermore, V"{qc) = V"{0) = u^, and 
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The evaluation of the determinant is not terribly difficult (the eigenvalues can be easily 
found; their product can be found in a table of mathematical identities); the result, for large 
(3, is 

From (1^), we can extract the ground state energy since, for large /3, Ke ~ exp —Eq(3/Tl. We 
find Eq = huj/2\iY> to corrections of order Xh^ . We have discovered an incredibly complicated 
way of calculating the ground state energy of the harmonic oscillator! 



9.3 Instantons in the Double Well Potential 

Let us now study a problem of much greater interest: the double well. We will see that 
configurations known as "instantons" make a non-perturbative correction to the energies. 
We wish to evaluate the PI 

Ke= / Pge-^^, 

where 

SE = jdr + ^{q' - a'f) . 

for (3 —>■ oo. As explained above, the PI is dominated by minima of Se, i-e., by classical 
solutions. The classical equation corresponds to a particle moving in the potential —V{q) 
(Figure 18); the "energy" E = ^q^ — V{q) is conserved. 

-V(q) 




Figure 18: Inverted double-well potential. 

Let us examine classical solutions, taking the boundary values q, q' of the classical solution 
corresponding to the maxima of —V, ±a. In the limit (3 — > oo, these will be solutions of zero 
"energy" , since as r — ±cx) both the kinetic and potential "energy" vanish. 

First, if q = q' = a (an identical argument applies if g = g' = —a), the obvious classical 
solution is g(r) = a; a quadratic approximation about this constant solution would be 
identical to the single-well case discussed above. 

But what if g = —a and q' = a (or vice- versa)? Then the obvious classical solution 
corresponds to the particle initially sitting atop the maximum of —V at —a, rolling towards 
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q — after a very long (infinite, in the limit /3 — > oo) time, and ending up at rest at the 
other maximum of — y as r — > oo (Figure 19). 

q(x) 

a 



Figure 19: Instanton in the double-well potential. 

We can get the analytical form of this solution: setting E ^ 0, we have 

There are a family of solutions interpolating between —a and a: 

q{T) = atanh ^(r - Tq), (58) 



where ui = and where tq is an integration constant which corresponds to the time 

at which the solution crosses q = 0. 

This solution is much like a topological soliton in field theory, except that it is localized in 
time rather than in space. One could argue that the solution doesn't appear to be locahzed: 
q goes to different values as r — > ±oo. But these are just different, but physically equivalent, 
ground states, so we can say that the instanton is a configuration which interpolates between 
two ground states; the system is in a ground state except for a brief time - an "instant" . 
For this reason, the solution is known as an instanton. 

I called the two solutions ^(t) = a and q{T) — atanh |-(r — tq) the obvious classical 
solutions because there are an infinite number of approximate classical solutions which are 
potentially important in the PI. Since the instanton is localized in time, and since the total 
time interval f3 is very large (in particular, much larger than the instanton width), a series 
of widely-separated instantons and anti-instantons (configurations interpolating between +a 
and —a) is also a solution, up to exponentially small interactions between neighbouring 
instantons and anti-instantons. Such a configuration is shown in Figure 20, where the hor- 
izontal scale has been determined by the duration of imaginary time /3; on this scale the 
instanton and anti-instanton appear as step functions. 

It is clear than an instanton must be followed by an anti-instanton, and that if the 
asymptotic values of the position are +a and +a the classical solution must contain anti- 
instanton-instanton pairs whereas if they are —a and +a we need an extra instanton at the 
beginning. 



9 INSTANTONS 



49 



Let us choose first limiting values g(— /5/2) = g(+/3/2) = +a. Then we are interested in 



In the spirit of (^4|), in the steepest-descent approximation Ke is equal to the sum of Pis 
evaluated about all classical solutions. The classical solutions are: gc(T) = a; Qc =anti- 
instanton-instanton= AI; = AIAI; etc., where the positions of the As and Is are not 
determined, and must be integrated over. Schematically, we may write 



where the superscript denotes the total number of Is or As. Let us discuss the first couple 
of contributions in some detail. 

qc = a: This case is essentially equivalent to the single-well case discussed above, and we get 



where to = is the frequency of small oscillations about the minimum of V. 

qc = AI: This case is rather more interesting (that is to say, complicated!). Let us suppose 
that the classical solution around which we expand consists of an anti-instanton at time ri 
and an instanton at T2 (Figure 21); clearly T2> Ti. 
Then we can write q = qc + y, and 



We can evaluate SE[qc]'- it is twice the action of a single instanton (assuming the I and A are 
sufficiently far apart that any interaction is negligible): 5'£;[gc] = 25'^'^'. The one- instanton 
action S'^''* is 




Ke = kI + kI + kI; + ---, 



(59) 




SeIq] = SE[qc] + S^'^'^ly]. 




q(x) 



-p/2 



p/2 ^ 



-a -- 



Figure 20: Multi-instanton configuration. 
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With the instanton profile given by (|58|) , the result is 



To evaluate the PI with the action S'^;"'^'^ [|/] , let us divide the imaginary time interval into 
two semi-infinite regions / and //, where the boundary between the two regions is between 
and well away from the A and the I (Figure 22). 

Then we can write 



K 



i?- ye 



1 1 +11 



(60) 



Here the first factor represents integration over the positions of the A and / (remember that 
the A must be to the left of the /!). The quadratic action can be written 



riquad riquad , riquad 

— I 



III 



where 5'^"'^'^^ is the quadratic action in the presence of an anti-instanton and 5'^"'^*^^-^ is that 
in the presence of an instanton. 

Then the PI separates into two factors: 



Ji+ii Ji Jii 



(61) 



where there is an implied integration over the intermediate position at the boundary of the 
two regions. The quadratic no-instanton PI also separates into two factors: 



Vye 



Vye~ 



■-iquadjO 



X / Vye'^E 
II 



quad.O 



./h 



(62) 



where the superscript "0" denotes that this is the PI about a no-instanton (constant) back- 
ground. We can combine ( |6TD and (^) to give: 



\,Vye-^T''\l^ ^^^Vye-^T^'^ul^' 



iquad 



Vye-^T'I"^ 



Vye 



(63) 



q(T) 
a.. 



Figure 21: Anti-instanton-instanton. 
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But 



and similarly for the last factor in (^), so we obtain 



(64) 



/+// 



V Tin 



where R is the ratio of the PI in the presence and absence of an instanton (or, equivalently, 
anti-instanton) given in (|64|) . Substituting this into (|60|), 



^2 ^ e-25rV^. 

^ V TT^l 2 



A similar argument gives 



^4 ^ g-45r A J^g-/3W2^4A^ 



4!' 



and so on for subsequent terms in the expansion 
Summing these contributions, we get 



1 + 



2! 



+ 



4! 



+ 



= y^e-^-/2^osh(/5i?e-^i^""/^ 

Now we must recall why we're calculating this object in the first place. The propagator can 
be written as in (^): 



Region / 



I 



Region U 



Figure 22: Division of imaginary time into two regions, one containing the anti-instanton, 
the other containing the instanton. 
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By comparing these two expressions we see that the lowest two energies are 

— - hRe-^'i''/^ and — + hRe-^'i''/\ 
2 2 

So the energy sphtting is given by 



AE = 2;^i^e-^-"/^ 



(65) 



AE is clearly non-perturbative: it cannot be expanded as a power series in h (or, equivalently, 
in A). 

In principle, we should calculate the ratio 

„ (instanton background PI) . r , 

R = — r- ~ ratio oi determinants, 

(constant background PI) 

but I don't know how to compute it other than by doing a very arduous, technical calculation; 
luckily, time will not permit it. The interested reader can consult the book by Sakita for a 
discussion of this calculation. 

As a final note, we have calculated the PI with q — q' — a; a good exercise is to do the 
analogous calculation for q — —a, q' = a. 

9.4 Instantons in a Periodic Potential 

Consider a particle moving in a one-dimensional periodic potential (Figure 23). 

V(q) 



q 

Figure 23: One-dimensional periodic potential. 

With two minima, as we have seen, instantons enable us to calculate the energy splitting 
between the lowest-energy states of even and odd parity. In a periodic potential, we will see 
that a continuum of energies arise. 

Let us label the classical minima of V by an integer, j. Clearly this model will have 
solutions analogous to the instantons above, going from any minimum of V to the adjacent 
minimum. We define an instanton as the classical solution going from any j to j + 1, and 
an anti-instanton as that going from j to j — 1. Then the Euclidean PI to go from j = to 
j = 0, for instance, can be computed in a manner similar to the calculation of the previous 
section. This time any number and any order of instantons and anti-instantons are possible, 
subject to the constraint that ni — ua- 
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A calculation similar to that of the previous section results in the following expression 
for the propagator: 



rrr 001/ 
V irh ^0 n!2 





nh n\ n'\ Jo 2n 

Jo 27rf-n n\ 




g 

n=U n'=0 



= yje-*"/^ r ^ exp {2/?i;e-«"/'' cos »} . 



The derivation of this is a worthwhile exercise. From it, we can read the energies: 

E{9)^^-2hRe-''^^''/^cos9. 

The second factor is an expression of the well-known result that the degeneracy is broken 

nonperturbativcly; the energies form a continuum, depending on the value of 6. 

Another application of instantons in quantum mechanics is the phenomenon of tunneling 
(barrier penetration). The instanton method can be used to calculate the lifetime of a 
metastable state in a potential of the form depicted in Figure 15. We will not discuss this 
application. 

Instantons also appear in (and are by far most useful in) field theory. In certain field 
theories the space of finite-Euclidean-action configurations separates into distinct topological 
classes. An instanton is a nontrivial configuration of this type. The necessary topological 
requirements for this to occur are not hard to satisfy, and the list of theories that have 
instantons includes the Abehan Higgs model in 1-1-1 dimensions, the 0(3) nonlinear a- 
modcl in 1+1 dimensions, the Skyrme model in 2-1-1 dimensions, and (most significantly) 
QCD. Instantons give rise to a host of interesting phenomena depending on the model, 
including confinement (not in QCD though!), 6'-vacua, a solution of the U(l) problem in 
strong interactions, and the decay of a metastable vacuum. Unfortunately time does not 
permit discussion of these fascinating phenomena. 



10 Summary and Gross Omissions 



In this set of lectures the subject of path integrals has been covered starting from scratch, 
emphasizing explicit calculations in quantum mechanics. This emphasis has its price: I have 
not had time to cover several things I would have hked to discuss. My hope is that having 
been subjected to calculations in gory detail for the most part in the relatively familiar 
context of quantum mechanics, you will be able to study more complicated and interesting 
applications on your own. 

Here is a list of some of the important aspects and applications of this subject which I 
didn't have time to discuss: 

1. Fermi fields and Grassmann functional integration; 

2. Gauge theories (gauge fixing and ghosts arise in a particularly elegant way); 

3. Feynman's variational method and application to the polaron (electron moving in a 
crystal environment); 

4. Derivation of the Landau-Ginsburg theory, including application to superconductivity; 

5. Instantons in field theory; 

6. Critical phenomena. 

I hope that, in spite of these unforgivable omissions, these lectures have been worthwhile. 
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